Saddle-node bifurcation cascade in optically injected lasers 
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Self-similar structures converging to accumulation points are observed in optically injected lasers. 
These structures are associated to saddle-node bifurcations. We found that these phenomena can be 
properly explained in the framework of saddle-node bifurcation cascade. Apparently not connected 
phenomena in these lasers are shown to be related within this framework. 
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Lasers with optical reinjection have played an impor- 
tant role in recent years, from three different points of 
view, namely theoretical, numerical and experimental. 
Recently, the presence of self-similar arrangements of pe- 
riodic orbits has been reported in this class of lasers 
[J, 0, 0, 0, 03]. The phenomena described in the liter- 
ature are associated to saddle-node (S-N) bifurcations 
and to periodic orbits accumulation points. Furthermore, 
smaller regions in parameter space have been reported 
showing the same behavior (see §3.2 in [li]). Also the 
presence of accumulation points of accumulation points 
has been described [2j]. 

A new mathematical framework is needed to organize 
most or all of these results, or rephrasing Q, to answer 
the question of what is the precise structure of the laser 
chaotic phases. In this letter we will show this framework: 
the saddle-node bifurcation cascade @, 0|. Formerly, we 
will show how the diverse phenomena described by re- 
searchers fit in this framework in a qualitative way. Later, 
after introducing the mathematical model of laser, we will 
show it in a quantitative way. This new framework would 
make the analysis of the problem easier, differentiating 
between basic processes and secondary effects. In fact, 
this model explains the following phenomena: 

(i) Self-similar structures. The self-similar structures 
above mentioned (associated to S-N bifurcations, with a 
fractal structure and converging to accumulation points 
that, on the other hand, converge to accumulation 
points), have the same structure as that found in the 
S-N bifurcation cascade, which was fully characterized in 

(ii) Accumulation of infinite cascades. A S-N bifurca- 
tion cascade is a sequence of S-N bifurcations in which the 
number of fixed points showing this kind of bifurcation 
is duplicated. If the whole canonical window is consid- 
ered, then infinitely many S-N bifurcation cascades exist, 
formed by sequences of points with periodicities q, q ■ 2, 
q ■ 2 2 , . . . , q ■ 2", . . . , for any n € N and q odd, called ba- 
sic period (see definition in §3.2 of [(|). These sequences 
converge, as n grows, to the same point (in a way resem- 
bling that of the accumulation horizons reported in 0). 
This point is known as the canonical window Myrberg- 
Feigenbaum point. Each of the q ■ 2™-periodic S-N orbits 



gives birth to a <7-2™-periodic window inside which all the 
process is repeated again as it was in the canonical win- 
dow, hence giving a fractal structure. Hence each of these 
g-2 n ~ windows has its own Myrberg-Feigenbaum accumu- 
lation point, and all these accumulation points converge 
toward the canonical window Myrberg-Feigenbaum accu- 
mulation point. In the injected laser literature this is re- 
ported by Bonatto and Gallas [3| as "boundaries formed 
by the accumulation of infinite cascades of self similar 
islands of periodic solutions of ever-increasing period" 
and "... these chaotic phases contain both single accu- 
mulations as well as accumulations of accumulations". 
The kind of phenomena described by Bonatto and Gallas 
matches perfectly what is called "attractor of attractors" 
in jf|, as we will see later in more detail. 

(in) Unnested period- doubling structures. In 0], 
Wieczorek, Kranskopf and Lenstra describe unnested 
period-doubling structures in injected lasers. They say: 
"the next bifurcation is the saddle-node of the periodic or- 
bits SL 2 in which two period-2 orbits, one attracting and 
one saddle are created". It is likely that they are describ- 
ing the two first elements of a S-N bifurcation cascade, 
that is, an unnested island of doubled period seems to 
be one element in a saddle node bifurcation cascade. In 
other words, the S-N bifurcation gives birth to a window 
that appears as an unnested island. 

(v) Intermittency. Wieczorek, Kranskopf and 
Lenstra also report chaotic behavior close to the S-N 
bifurcation (jj. They are observing one single element 
of the intermittency cascade associated to the S-N 
bifurcation cascade In summary, before the S-N 

bifurcation, there is intermittency, and after it the 
system enters a periodic window where the initial 
periodic orbit undergoes a period-doubling cascade 
leading to chaos. We will show this intermittency below. 
Within the theoretical framework just described, the 
S-N bifurcation is the main element and the secondary is 
the intermittency which is, in fact, naturally associated 
to it. 

It is enough to find a S-N bifurcation cascade in rein- 
jected lasers to apply the theoretical framework just de- 
scribed to these systems. In this letter, we find S-N bifur- 
cation cascades in the same reinjected laser model used 
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by [H, 0, S 3 j providing the mathematical basis to ex- 
plain the phenomena described above. We have chosen 
this laser model because periodic orbits associated to S- 
N bifurcations have been found in it which agree with 
experimental results. Furthermore, the model is quite 
general and representative of a broad number of class-B 
lasers (solid state and CO2 lasers), as pointed out in []]]. 
In this way, we expect to ease the experimental research 
based in our theoretical developments. The equations of 
the model can be written in a dimensionless form as 

E = K+(-(l + ia)n - iu))E 

h = -2Tn - (1 + 2Bn)(\E\ 2 - 1) 

(see [H for details) ; the parameter values have been taken 
as a = 1.987, uj = 1.5, T = 0.035 and B = 0.015. 

Taking K as the control parameter, S-N bifurcation cas- 
cades are found (see figure [J) . In tables HI [ill and 1111} the 
values of K where S-N orbits of several periods appear 
are shown, as well as the Feigenbaum ratios, converg- 
ing to the Feigenbaum constant 5. They show the 3 • 2™ 
and 5 • 2™ cascades as well as the S-N bifurcation cascade 
within the period-5 window. This proves, in optically 
injected lasers: 

(a) The existence of cascades with different basic peri- 
ods (3 and 5). This explains the existence of unnested 
islands. 

(b) The same kind of behavior in smaller regions (S-N 
bifurcation cascade in the period-5 window). This ex- 
plains the existence of self-similar structures. 

(c) The underlying attractor of attractors originated 
by the presence of cascades inside each of the windows of 
the cascade. This explains the accumulation of accumu- 
lation points. 

Then, the phenomena reported by researchers are ex- 
plained within the S-N bifurcation cascade framework. 

Let us point out that all of the S-N bifurcation cas- 
cades have the same scaling as the Feigenbaum cascade 
0, 0| allowing this property for a technique to have them 
characterized by researchers (see tables [U HU and . 

Another point has to be highlighted, the g-periodic or- 
bit, for q odd, is located in the 1-chaotic band, the q ■ 2, 
lies in the 2-chaotic band, and so on. This allows us to 
explain the structures described by Bonatto and Gallas. 
For instance, orbits 3, 5, 7 (the last one not numbered 
but shown) in Fig. lb in 0, belong to the 1-chaotic or- 
bit. The elements located in the 2-chaotic band allow to 
explain the low periodic islands with periods 10, 14, 18, 
22 shown in Fig. 2b in Q. They are low periodic islands 
of periods 5-2, 7-2, 9-2 and 11-2 which are respectively as- 
sociated with the second elements of the S-N bifurcation 
cascade q-2 n , having basic periods q = 5, 7, 9 y 11. The 
periods of the islands just mentioned follow a Sharkovsky 
ordering [13], which claims that the low periodic island 



of period 6 (shown in figure lb in Q|) is in reality a 3 • 2- 
periodic island, that is, it belongs to the 2-chaotic band 
and represents the second element of the S-N bifurcation 
cascade with basic period q = 3. Closeups of these peri- 
odic islands can be seen in the biparametric plot in figure 
[3] (similar to figure lb in Q). If the biparametric struc- 
ture is visited along a straight line across the unnested 
periodic islands, then a plot similar to the one in figure Q] 
is obtained from which the S-N bifurcation cascades can 
be extracted. 

In this same framework, it is possible to explain the 
two kinds of convergence reported by researchers: 

(d) Simple accumulation. If for q ■ 2™, n is kept fixed 
and q is increased, the corresponding sequence repre- 
sented by the low periodic islands in the 2™-chaotic 
band converges to the Misiurewicz or band-merging point 
where two consecutive bands meet each other [t| . This is 
seen in Fig. 2a in [2] where n = 1 and q = 5, 7, 9, 11, . . . 
and the low periodic islands converge to an "accumula- 
tion horizon". 

(e) Accumulation of accumulations. If, for a fixed q, 
n is increased, the corresponding sequence converges to 
an attractor of attractors (see [6j). 

In fact, all the convergence points are attractors of at- 
tractors when observed in enough detail, because within 
any window there exists an infinity of attractors of at- 
tractors. However the software or experimental setups 
employed cannot resolve the simple accumulation and 
therefore the accumulation of accumulations cannot al- 
ways be observed. 

Finally, there must be an intermittency before any S- 
N bifurcation. In particular, 0| report an intermitency 
behavior. But, in general, there must be intermittency 
surounding every low q • 2™-periodic periodic island de- 
scribed above. These would be responsible of the appear- 
ance of intermittency cascades de intermitencias behav- 
ing in a similar way as S-N bifurcation cascade do. We 
show one of these intermittencies in figure [H 

According to the exposed in the previous paragraph, 
to show up the existence of a S-N bifurcation cascade 
in a dynamical system the chaotic band containing ev- 
ery q ■ 2"-periodic orbit (for a fixed q basic period) of 
the cascade has to be perfectly identified. Otherwise, 
that is, if the succesive windows that appear when the 
control parameter is varied are registered in order, then 
the S-N bifurcation cascade will not be readily obtained, 
but a mixture of unrelated periodic windows, that do not 
follow the Feigenbaum scaling any more. The S-N bifur- 
cation cascades are subsequences of this sequence that 
have to be chosen carefully according to the (symbolic) 
rules given in 0|. 

In summary, all kind of phenomena (self-similar struc- 
tures, accumulation points, unnested period-doubling 
structures, intermittency and so on) reported by re- 
searchers in optically injected lasers can be explained 
in the theoretical framework of saddle-node bifurcation 
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FIG. 1: Fragment of the period-5 window. Several periodic 
windows are seen within it. The origin of each of these pe- 
riodic windows (its right bound) coincides with the birth of 
saddle-node orbits of the corresponding saddle-node bifurca- 
tion cascade. 
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FIG. 2: Intermittency behavior close to the onset of the 5- 
periodic window, for a value of k — 0.938818. 



cascade. Furthermore, new unreported phenomena (in- 
termittecy cascade) have been predicted. 



FIG. 3: Low periodic islands of different periods (3-2 and 5-2) 
born at S-N bifurcations and located in the 2-chaotic band. 
These S-N bifurcations belongs to different S-N bifurcation 
cascades (basic periods 3 and 5 respectively). The points in 
this plot have been computed by period-counting. 



period 




(«„_! - Kn-l)/(K n - Kn-l) 


3x2 


0.984222 




3 x 2 2 


0.994836 




3 x 2 3 


0.997134 


4.619 


3 x 2 4 


0.997621 


4.719 


3 x 2 5 


0.997726 


4.638 



TABLE I: Values of k corresponding to S-N bifurcation cas- 
cade of period 3 in the canonical window. 



period 




(Kn-l - Kn-a)/(«n ~ «n-l) 


5 


0.938815 




5x2 


0.989430 




5 x 2 2 


0.995914 


7.806 


5 x 2 3 


0.997363 


4.475 


5 x 2 4 


0.997671 


4.705 



TABLE II: Values of k corresponding to S-N bifurcation cas- 
cade of period 5 in the canonical window. 



period 




\f^n—l K>n — 2)/(^n K>n — l) 


(5 x 3) 


0.937958 




(5 x 3) x 2 


0.9380750 




(5 x 3) x 2 2 


0.9380993 


4.815 


(5 x 3) x 2 3 


0.93810450 


4.673 



TABLE HI: Values of K corresponding to S-N bifurcation cas- 
cade of period 3 inside the period-5 window. 
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